Drift Lagrangian for relativistic particle in intense laser field 
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The Lagrangian and Hamiltonian functions describing average motion of a relativistic particle 
under the action of a slightly inhomogeneous intense laser field are obtained. In weak low-frequency 
background fields, such a particle on average drifts with an effective, relativistically invariant mass, 
which depends on the laser intensity. The essence of the proposed ponderomotive formulation is 
presented in a physically intuitive and mathematically simple form, yet represents a powerful tool 
for studying various nonlinear phenomena caused by interaction of currently available smooth ultra- 
intense laser pulses with plasmas. 

PACS numbers: 52.35. Mw, 52.20.Dq, 52.27.Ny, 45.20.Jj 



43 

Oh! 



CZ3 

a _ 
11; 

o ■ 

• I— I . 

>v 
43 ' 



> ■ 

Os ■ 

o . 

o ■ 
cn : 
o . 

o : 

43 : 

Or 



X 



The latest advances in high-power laser technology 
have resulted in the development of laser systems capable 
of delivering superstrong electromagnetic pulses, which 
can be focused to intensities as high as 10 21 W/cm 2 p|, 
with even more powerful systems coming up in the near 
future [2j. The currently obtainable laser helds can ac- 
celerate electrons up to ultra-relativistic oscillatory ve- 
locities, previously unachievable in experiments on laser- 
matter interaction. This revolutionary progress is now 
giving a new life to theoretical studies on particle be- 
havior under the action of intense electromagnetic radia- 
tion. The conventional models describing various nonlin- 
ear phenomena in plasmas illuminated by high-frequency 
radiation nowadays need to be revised, as new, ultra- 
powerful laser systems are becoming available for lab- 
oratory experiments. To explain the already observed 
phenomena and predict new effects taken place under 
the action of intense laser drive, adequate description of 
single-particle motion under relativistically strong radia- 
tion must be developed first. 

Currently, particle motion is well-understood when the 
only forces present are those from the wave of uniform in- 
tensity 0. However, to study the guiding center dynam- 
ics in inhomogeneous laser radiation or drifts determined 
by the presence of low-frequency background fields, ad- 
ditional analysis is needed. Expanding the conventional 
understanding to this area would result in a substan- 
tial progress in studying a number of plasma physics 
problems, such as, e. g., Coulomb collisions and ener- 
getic particle production in strong laser fields. Moreover, 
the hydrodynamics and the electrodynamics of laser- 
illuminated plasmas would readily be available for gen- 
eral revision. 

Often, the dynamics of a particle moving in a high- 
frequency field is described in terms of the ponderomotive 
approach. In nonrelativistic ponderomotive description, 
the effect of high-frequency electromagnetic forces on a 
particle is replaced by particle interaction with an aver- 
age potential, linear on the intensity of laser radiation 
0| . When ultra- intense lasers are employed, this conven- 
tional description needs to be generalized to relativistic 



motion. Contrary to the degenerate case of a circularly 
polarized field, in which high-frequency variations of the 
relativistic mass can be neglected [B| , the problem of par- 
ticle motion in the case of a linear or elliptic polarization 
represents a certain challenge, but still can be studied 
analytically. 

To describe the drift particle dynamics in such fields, 
recently, multiple studies were performed. Under various 
approximations, it was shown that the oscillating par- 
ticle guiding center drifts in a smooth laser field with 
an effective mass, which depends on the electromagnetic 
field intensity [E E IE • (In strongly nonuniform laser 
fields though, the particle dynamics is more complicated 
[E inj.) The Hamiltonian treatment of the relativis- 
tic drift under intense laser drive has been proposed in 
Ref. though the problem of interaction with low- 

frequency background fields has not been studied. The 
first steps towards developing the general formalism with 
such interaction were made in Ref. ||. However, only 
smooth (compared to the amplitude of oscillations) low- 
frequency background fields were taken into considera- 
tion and the relativistic drift motion equations were in- 
duced without proper justification. 

These shortcomings are overcome in our paper, which 
major emphasis is twofold. First, we propose a general, 
fully relativistic Lagrangian formulation of ponderomo- 
tive description of particle motion under the action of a 
quasi-monochromatic slightly inhomogeneous laser field. 
The proposed approach is physically intuitive yet more 
systematic and simple in comparison with those discussed 
previously. After natural generalization, it allows includ- 
ing particle interaction with weak background forces, ad- 
ditional to those from the laser field. The discussion on 
that aspect of the guiding center motion constitutes the 
second emphasis of our work. We show the effective mass 
concept to be applicable to ponderomotive description of 
relativistic particle motion in low-frequency background 
fields, including even ones of small spatial scale compared 
to the amplitude of oscillations. In the end, we discuss 
the most promising applications of the proposed formu- 
lation and summarize our main ideas. To start, consider 
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particle motion under the action of a plane laser wave 
propagating in vacuum, with vector potential given by 

A(r,t) = (mc 2 /e) 0(77), (1) 

where 77 = tot — k ■ r stands for the phase of the wave, 
u) is the wave frequency, and k = z°lu/c represents 
the wavevector. The polarization of the wave will be 
assumed fixed though arbitrary. The magnitude of a, 
a = eE/mao (where E is the laser electric field), can be 
understood as the ratio of the momentum imparted by 
the wave field in a single oscillation to mc, meaning that 
relativistic effects become important at a > 1. (For the 
wavelength A = 2ttc/uj = 1 /im, the intensity correspond- 
ing to a ~ 1 for electrons is about 10 18 W/cm 2 .) 

In a certain, unique, frame of reference, in such a 
field the particle undergoes stationary oscillatory "figure- 
eight" motion in a linearly polarized wave or circular mo- 
tion in a wave with circular polarization Averaging 
over the oscillations, one comes to the concept of the 
guiding center motion, which we study below. First, let 
us consider the variational principle that states the min- 
imum value of the action 

S = f 2 Ldt, (2) 

where L is the Lagrangian function of the particle mo- 
tion to be realized on the true trajectory. On time scales 
ti — t\ large compared to the oscillation period, the ma- 
jor contribution to the action S (linear on £2 — t\) is 
provided by the time-averaged part of the Lagrangian, 
(L) , while the contribution of the oscillatory Lagrangian 
into the integral @ remains small. Thus, the action S 
is approximately given by S = f t 2 {L)dt, from where it 
follows that (L) can be treated as the Lagrangian of the 
average, guiding center motion. 

To obtain the form of the drift Lagrangian Lq = (L) , 
let us consider the latter in the frame of reference where 
the guiding center rests. In the new frame, the guid- 
ing center Lagrangian Co can be nothing but a constant, 
which we put in the form 

C = -m cS c 2 , (3) 

in analogy with the Lagrangian of a true particle with 
zero velocity. The formally introduced quantity m e a 
playing a role of a new, effective mass is yet to be de- 
fined. The action J2J) is relativistically invariant and can 
be written as S = Cq dr, where the time r represents 
the proper time of the guiding center. Since dr is in- 
variant by definition (and thus, so is the Lagrangian Cq), 
the quantity m e ff must also be relativistically invariant. 
Using 



dr = dtJl- v 2 /c 2 , Co dr = L dt, (4) 



guiding center motion 



L = -m ffc 2 y I - vl/c 2 , (5) 

which formally coincides with the Lagrangian of a rela- 
tivistic particle with mass m e e moving with velocity Vo- 
Since the original frame was chosen arbitrarily, the above 
expression represents the general form of Lq, where m e g 
is left to be expressed in terms of the parameters of the 
laser field. 

Let us calculate Lo in a laboratory frame of reference 
where the particle has a nonzero average velocity v - In- 
structive by itself, the derivation to follow will also pro- 
vide us with a number of useful relations connecting the 
parameters of the particle drift and those related to the 
actual motion. To proceed, consider the Lagrangian of 
the particle true motion given by 



L = -mc^l-^ + - c (vA( V )), (6) 

which is a known periodic function of the phase 77 rather 
than time t. Thus, in order to average L over time, one 
needs to derive a relation connecting time averaged and 
phase averaged quantities. For an arbitrary quantity /, 
its time and phase averaged values given by 

f + A )) + 27T 

</> = £ / fdt', J=± J fdrf, (7) 
t 77 

where the limits of integration over the phase correspond 
to the limits of integration over the time (i. e. 77 = rj(t)), 
and the time interval A is defined as one on which the 
total phase change equals 2ir: 

rf+2-K 
V 

The time interval A coincides with the wave period 2ir/u 
only if particle motion is nonrelativistic. However, gen- 
erally, the phase time-derivative is given by 

dt =UJ { 1 ~-)= UJ { 7 J' (9) 

where 7 = (1 — v 2 /c 2 ) is the normalized relativistic 
energy [E — mjc 2 ) and p — mjv it the particle kinetic 
momentum. 

Since the original Lagrangian depends on 77 (that is, on 
the combination z — ct, rather than z and t separately), 
there exists an invariant of motion given by 

u = 7 — p z /mc — const. (10) 

Substituting the above expressions into Eq. Q , one gets 



where v is the velocity of the guiding center in the orig- 
inal frame of reference, one gets the Lagrangian of the 



(/) = 7//T 



(11) 
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Note that the obtained formula is valid only in case when 
the electromagnetic wave Q is propagating in vacuum. 
If the refraction index of the medium differs from unit, 
Eqs. © and i|10[) need to be modified, and the relation 
between time- and phase-averaged quantities becomes 
more complicated jl 1| . 

From Eqs. (HUJ and ITTJ) , it follows that 

7= y/l + (p/mc) 2 +tf, 7 = 70^^+^, (12) 
where 70 = (1 — Vq/c 2 )^ 1 / 2 , and 

v = (v) = p/rrvy (13) 

is the drift velocity of the particle (compare with the 
inexact expression given in Ref. Q). Thus, Lq can be 
put in the form (jSJ) with m e ff given by 



m ofr = m\j 1 + e 2 A 2 /m 2 c 4 . (14) 

The guiding-center Lagrangian JSJ with the expression 
(|14fl for the effective mass was also obtained in Ref. 9] by 
a somewhat similar yet a complicated and not a straight- 
forward procedure. In the cited work, Eq. (|14|) was sup- 
posed valid only in the frame of reference where Vq = 0. 
In fact, as shown above, it remains applicable for arbi- 
trary Vq, and, more than that, the actual value of m e g 
must be relativistically invariant. To express the effec- 
tive mass in the invariant form, let us notice that, in the 
laboratory frame where we chose the electric potential 
4> = (see Eq. ©), VA 2 coincides with the norm of 
the 4- vector potential V ' A a A a , A a = (<f>, A). The latter 
is Lorentz- invariant 3] , and remains such after being av- 
eraged over relativistically invariant phase rj. Thus, the 
expression for m e ff, invariant to relativistic transforma- 
tions, can be put in the following form: 



m cS =m\ l+ — {A a A a ). (15) 

Eq. (|15fl was also given in Refs. 8] where the average 
particle motion was studied otherwise. 

Reverting to the formula for the drift Lagrangian (JSJ) 
with the effective mass given by l|14[l . the canonical mo- 
mentum of the guiding center motion Pq equals the 
phase-averaged kinetic momentum p = m e s"foVQ, and 
thus the Hamiltonian function of the guiding center mo- 
tion can be put in the form 

H = ^ml^ + P 2 c 2 . (16) 

Here m c ff may smoothly depend on the guiding center 
location Rq and time t if the wave envelope is slightly 
nonuniform or time-dependent. Precisely, that means 
that the laser intensity "seen" by the particle changes 
insignificantly on one period of particle oscillations, so 
that the averaging still makes sense, i .e. 

;>r^, T>A, Z/u >A, (17) 



where I and T are the spatial and the temporal scales of 
the wave envelope (for detailed analysis, see Refs. [rllTol 

E2). 

An alternative derivation of Eq. (|16fl c & n be found in 
Ref. , where a sequence of canonical transformations 
of the original motion equations was shown to lead to 
a similar result. In the present paper, we showed this 
tedious procedure to be unnecessary for obtaining the 
expression for the drift Hamiltonian. Comparing to the 
cited work, the distinguishing advantage of the formula- 
tion proposed in the present paper is that, because of 
its apparent mathematical simplicity, this formulation 
allows easy generalization of the drift Lagrangian and 
Hamiltonian formalism on the case when the oscillating 
particles undergo weak acceleration by large-scale low- 
frequency forces satisfying l|17|) . Interaction with these 
forces enters the expression for Lq additively and, what 
is most important, can still be considered in the frame- 
work of the effective mass concept. 

To show this, consider an oscillating relativistic parti- 
cle interacting with a field governed by the 4-vector po- 
tential A^ = (</>bg, ^4.bg), where the subindex "bg" stands 
for a background field, additional to the one of the laser 
wave. Assume that the field is weak: 

eE hg /j m cS c < u, eB hg /j m cf fC < uj, (18) 

where .Ebg and _Bb g are the corresponding electric and 
magnetic fields. In this case, the background fields do 
not impact the oscillatory motion significantly. Thus, 
averaging of the kinetic term mc 2 /7 in the Lagrangian 
leads to the same expression as in Eq. (JSJl with m g- given 
by Eq. I|15|l . In the zeroth-order approximation with 
respect to the small parameters (|1S|> . the average part 
of the Lagrangian corresponding to particle interaction 
with the background field can be expressed in terms of 
the quantity Aq = (</>o, Aq) given by 

A * = (A% g (R +r^). (19) 

The time-averaging procedure is invariant with respect 
to changing the drift frame of reference, i. e. does not 
alter the Lorentz transformation properties of the quan- 
tity being averaged. Thus, Aq represents a true 4-vector 
and can be considered as a new, effective electromag- 
netic field. In terms of this field's potentials, the drift 
Lagrangian can be put in the following form: 



Lo = -m cS c 2 ]Jl--^ + --(vo-Ao)~e(j)Q. (20) 

In certain applications, it is of interest to consider par- 
ticle interaction with background fields having spatial 
scale /bg S r ~- If the drift velocity is small, so that the 
drift displacement on a single period uoA is small com- 
pared to ?bg) the ponderomotive description still can be 
applied. However, in this case the difference between the 
time-averaged potential Aq and the true potential A" g 
taken at the location of the guiding center Rq, is crucial. 
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For example, this situation is realized at Coulomb scat- 
tering in intense laser fields when exceeds the radius of 
effective interaction Note that, as follows from the 
above analysis, the characteristic amplitude of the effec- 
tive potential remains unchanged as one generalizes the 
expression for </>o to the case of rclativistic particle mo- 
tion. In this case, the only difference in calculating <j>o is 
provided through the change in the oscillatory trajectory 
r~(t) to be averaged over. 

In the context of the Coulomb scattering problem, the 
considered Lagrangian approach represents a unique tool 
for studying ponderomotive and even stochastic behavior 
of particles being scattered. This problem deserves de- 
tailed consideration and will be discussed in future works, 
though, briefly, the extension of the proposed formula- 
tion can be explained as follows. Stochastic behavior 
of a dynamical system with periodic coefficients is often 
convenient to describe in terms of mapping of the dynam- 
ical trajectory onto a subspace of the system phase space 
(for review, see Ref. [T^l. For the Hamiltonian mapping 
(R ,P ) — > (Ro,Po) connecting the particle locations 
and momenta before and after the time interval equal 
to the period of the laser field, the generating function is 
given by the action with t± =t and t 2 —t + 2n/uj [l4|. 
Since the obtained drift Lagrangian is approximately pro- 
portional to S, it can readily be used for constructing 
the actual form of this mapping. As will be shown in our 
future publications, when studying the statistical prop- 
erties of particle stochastic dynamics (rather than sin- 
gle particle motion) by means of such a mapping, the 
conditions l|17|) can be significantly relaxed, which sub- 
stantially broadens the proposed Lagrangian approach 
applicability. That also allows a significant progress in 
studying the problem of energetic particle production in 
strong laser fields [l5| . 

Since, in the case of relativistic drift, depends on 
Vo, the expression for the canonical momentum Po — 
dLo/dvo the drift motion equations become complicated. 
However, in two special cases of interest, those can be 
simplified. In a large-scale background field satisfying 
the conditions (|I7J1 , locally, A" g can be treated as a linear 
function of r. Therefore, the velocity-dependent part 
averages out when calculating the potential Aq , and one 
gets A% b 

equals P = m cS ~/ v 
function is given by 



A^ £ . Thus, the drift canonical momentum 
(e/c) -4.bg, and the Hamiltonian 



H = Jm 2 eff c i +[Po--A hg ) c 2 + e<p hg , (21) 



where the potentials are assumed to be slow functions of 
Ro and t. The guiding center motion equations can be 
put in the covariant form 

dRo _ Po dp% _ e aP _ 2 dm cf f 

dr ~ m eff ' dr ~ c b « Up dR* ' [ ' 

where Rq = (ct, Rq) is the 4-coordinate of the guid- 
ing center, Pq — (£o/c, m e fT7oVo) is the drift kinetic 



4-momentum, £ = m e ff7 c 2 is the energy of the guid- 
ing center motion, is the electromagnetic field ten- 
sor corresponding to the potential A% Q, and U a = 
70 (c, Vo) is the guiding center 4- velocity. Covariant Eqs. 
(|22|l were also given in Ref. though no strict deriva- 
tion of those was proposed. Another expression for the 
relativistic ponderomotive force is given in Ref. [jj]. 

From Eq. I|2II) . it follows that, in low- frequency large- 
scale background field, the guiding center of a relativistic 
particle moving under the action of intense laser radiation 
behaves as a particle with the effective mass m e e drifting 
in the same background field. For example (and in co- 
incidence with the results obtained in Ref. 8]), in static 
magnetic field -Bbg, the guiding center undergoes Larmor 
motion with the cyclotron frequency u>b = e-Bbg/7ow e ffC. 
Conventional expression for the drift velocity in nonuni- 
form magnetic field Q also readily applies to the average 
motion if the particle true mass is replaced with the one 
given by Eq. Q15[l. 

In addition to the case of large-scale background fields, 
the guiding center motion equations can also be put in 
a simple, physically intuitive form in the case of nonrel- 
ativistic drift motion. Since the drift velocity enters the 
expression for Aq only through relativistic dependence of 
on Vq/c, then, in the case v <C c, OAq/Ovq can be 
neglected. In this case, the drift canonical momentum is 
given by Po = WcffVo + (e/c)Ao, and the Hamiltonian 
can be put in the form 

H = 7T— ( p o --A ) + m cff c 2 + e0 o , (23) 

where the effective mass m e ff and the potential energy 
ip e g = m c gc 2 + etfio may slowly depend on the guiding 
center location Rq and time t. Note that even in uni- 
form laser field, VV> e ff may differ significantly from eV0bg 
when the amplitude of particle oscillations exceeds the 
spatial scale of the background field I . The regime of 
slow drift motion superimposed on relativistic oscillations 
is the one, which is actually realized in many current ex- 
periments on intensive laser pulses interaction with rare 
plasmas. (Under rare plasmas we assume those having 
a refraction index close to unity, as assumed for all the 
results obtained in the present paper.) This fact makes 
the above analysis especially useful from the practical 
point of view, as it represent a simple tool for studying 
the actual experimental data. Finally, the well-known 
nonrelativistic ponderomotive potential can be read- 
ily derived from Eq. ll'-il) by keeping the correction to the 
effective mass, linear with respect to the wave intensity 
(see also Ref. 0)- 

In summary, we showed that, in weak low- frequency 
background fields, relativistic particle moving under the 
action of intense laser radiation drifts like a quasi-particlc 
with an effective mass, which depends on the intensity of 
the laser field. The intuitive expectation that, by the 
order of magnitude, the drift motion equations must co- 
incide with those without the laser field if the appropriate 
relativistic correction of particle mass is introduced, can 
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now be considered proven for various types of background 
fields. The proposed formulation can be useful for study- 
ing numerous phenomena resulting from intense laser- 
plasma interaction, such as, e. g., the energetic particle 
production and Coulomb scattering in strong laser fields. 
Moreover, the mathematical simplicity of the proposed 
approach allows easy generalization of the rare plasma 
hydrodynamics and electrodynamics to the case of plas- 



mas illuminated by ultra-intense laser radiation. Replac- 
ing the electron mass with the effective mass (|15f) . one 
can readily derive the generalized dispersion relations for 
various linear waves in plasmas, as well as revise the non- 
linear plasma dynamics. 
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